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Abstract 

We use the 2PI effective action of a relativistic scalar field theory to 
derive kinetic equations for a Bose-condensed system near the phase 
transition. We start from equations of motions derived within a jt- 
expansion at NLO. In taking the non-relativistic limit we obtain a 
generalized Gross-Pitaevskii equation for the condensate field. Within 
the Popov approximation we explicitly compute the collision term up 
to order g 2 using the Kadanoff-Baym formalism. For the sake of self- 
consistency we derive in the same way a Boltzmann equation for the 
non-condensate distribution function. The final results are in agree- 
ment to those previously obtained by Griffin, Nikuni and Zaremba 
and by Stoof. 
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1 Introduction 



Since the first experimental realization of a Bose condensed dilute gas in 
1995 there has been an intense effort to further understand this new state 
of matter both from the experimental and theoretical point of view. The 
most common theoretical ansatz for a wide range of questions is the famous 
Gross-Pitaevskii equation of motion (e.o.m.) JTj. This equation for the Bose 
macroscopic wavefunction $ gives a complete description of the condensate at 
T = 0. Near the phase transition one clearly has to take into account the ef- 
fects of the non-condensate atoms. One approach to this problem is to derive 
a useful kinetic equation that includes the collisions between the condensate 
and non-condensate part of the system. This was done by Griffin, Nikuni and 
Zaremba [2J starting from a suitably generalized Gross-Pitaevskii equation 
(see also the work by Imamovic-Tomasovic and Griffin [31 EJ E])- Further 
approaches to a quantum kinetic description of Bose condensed systems can, 
e.g., be found in El H • 

In this paper we apply a method coming from non-equilibrium relativistic 
quantum field theory. For a recent review of the so called two particle ir- 
reducible (2PI) effective action techniques see ^Oj- This powerful method 
has already been used in a variety of interesting issues like physics of the 
early universe and heavy ion collisions fTJ EI] , and recently to the non- 
equilibrium dynamics of Bose-Einstein condensates (Hj. We use the e.o.m. 
for a N-component scalar field obtained within a -^-expansion of the 2PI 
effective action to next to leading order (NLO) ^H]- We finally arrive at the 
same kinetic equation previously derived in E] ■ 

2 Field equation of motion 

Our starting point is the relativistic e.o.m. from the -^-expansion of the 2PI 
effective action at NLO. For a N-component real scalar field <p a with mass m 
and |y0 4 -interaction one has (see (27) in [TK] ) 

- (o + m 2 + A [$ c (x)$ c (z) + G cc (x, x)]j $ (x) = K a {x, x), (1) 
with (see (35) in [ISj ) 



K a (x,y) = -^<S>b(x)G ba (x,y) 



_ *g^ J dz ° J d 3 zG bc (x,z)G bc (x,z)K a (z,y), 



(2) 



2 



and the propagator G given by 



G ab (x,y) =< T(p a (x)(f) b (y) > - < (f> a (x) >< <j> b {y) > . (3) 

The mean field is denoted by <& a (x) =< <f> a (x) >. In (3) T means time or- 
dering along the Schwinger-Keldysh contour ^HJEl^lEII]- The z°-integral 
in (2) is taken along this contour as well (see fig.l). 

As we ultimately want to describe a Bose condensed gas we have to take 
N = 2, or a = 1,2. 

Next we perform the non-relativistic limit by the following replacements 

$ a — -U^* (4) 



G — > —G NR , (5) 
m 



n — -Ln™, (6) 



K a — 2^K» R , (7) 
A — ► 12m 2 g, (8) 

(□ + m 2 ) -^-2m(^o + — ), (9) 

2m 

n(x, y) = --Ga&O, y)G ab (x, y). (10) 

The dimensionful prefactors arise because of the different normalizations of 
relativistic and non-relativistic fields. Throughout the paper we set h = c = 
1. 

The non-relativistic e.o.m. now reads 



where we define 



id x0 + a - 1 p(x)$f R (x)+c(^)]) = 

g^ R (x)Gil R (x,x) + 2ig J^dz J dh U NR (x, z)K^ R (z, x). (11) 
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For a later comparison to other work on Bose condensed systems it is conve- 
nient to express (11) in terms of one complex field <£>. Thus we perform the 
following redefinition (we suppress from now on the 'NR' index) 



_>-!(¥* + *), 2 ^_L(**-*), (12) 

with 

^ = $ + ^, (13) 
and correspondingly for the expectation values 

< ^ >= $, <v£> = 0. (14) 

The interpretations of <£> as the condensate and \1/ as non-condensate complex 
fields are obvious. The components of the propagator matrix G a b are replaced 
by 

1 

G n (x,x) — > -(to* + 2n + m), (15) 
G 2 2(x,x) — > ~(rh* - 2h + to), (16) 

l 

G 12 (x,x) — >-(m*-m), (17) 

G 2 i(x,x) — G± 2 (x, x), (18) 
where we introduced the non-equilibrium non-condensate densities 

n( x ) =< ^*(x)^(x) >, m(x) =< ^(x)^(x) > . (19) 

The e.o.m. (11) reads in terms of the new fields 

id x o H g [n c (x) + 2n{x)\ J $(x) — gm{x)Q*{x) = 

2m J 

2ig [ dz° [ d 3 zU(x,z)K(z,x), (20) 

with 



c 



K(x,y) = -j=(K 1 (x,y)+iK 2 (x,y)) 

= g$(x) < TV*(x)V(y) > +g$*(x) < TV(x)V(y) > 

+ 2ig J dz° J d 3 zIl(x,z)K(z,y), (21) 
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and the non-equilibrium density in the condensate n c (x) = |($^(a:)-|-$|(a;)) = 
|$(x)| 2 . 

It is convenient to write the r.h.s. of (20) in terms of IT^ and K<. These are 
defined as EH 



n > (x,?/) = Tl(x,y) for x° > y° 



n < (x, y) = y) for x°<y°, 

and equivalently for K<. 
Using the relation 



(22) 
(23) 



J dz° J d 3 z U(x,z)K(z,x) = 

J dz° J d 3 z(U > (x,z)K < (z,x) -U < (x,z)K > (z,x)) (24) 
finally leads to the e.o.m. 



(id x o + ^ — g[n c (x) + 2n(x)])&(x) — gm(x)§* (x) = 
2ig J dz° J d 3 z{U > (x,z)K < (z,x) -U < (x,z)K > (z,x)), (25) 
together with (21) for K<. 

This equation can be seen as a generalization of the time dependent Gross- 
Pitaevskii equation. In particular, the RHS contains interactions between 
the condensate and the surrounding non-condensate cloud. 





ilmz 
















Rez° 



Figure 1: Schwinger-Keldysh contour in the complex 2°-plane 
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2.1 Derivation of the kinetic equation 

Our goal in this section is to derive a kinetic equation for $ in terms of 
the non-condensate distribution functions. This equation will have a clear 
physical interpretation. We have to make some approximations as the very 
concept of a distribution function f(k,x,t) cannot be strictly applied to a 
quantum mechanical problem. First of all, we will work within the Popov 
approximation [20J, i.e. we neglect the anomalous correlators < > and 
< >. In addition, we suppose the coupling constant g to be small as 

in current experiments the densities near the phase transition are such that 
the gas is in the weak coupling limit [.[}. 
One then obtains 

U^(x,y) ~ - < Of(x)V*(y) >^< V*(x)W(y) (26) 



K^(x, y) ~ g$(x) < ^*{x)^(y) (27) 

with < ^(x)^*(y) ><=< Sk*(y)iff(x) > and so on. We take only the first 
term of the iterative equation (21) for K< what leads to a collision term up 
to order g 2 . 

Following the well known approach of Kadanoff and Baym [21], we now 
assume that the correlation length of the system under consideration is suf- 
ficiently small and one can write 

< >= J ^eV^f&X), (28) 

< %)*%) >= J -^e**"*" (1 + f(p, X)) , (29) 

where s = x — y and X = are the relative and center of mass coordinates, 
respectively, while p — (u,p) denotes the non-condensate quasiparticle four- 
momentum. The dependence of the correlation functions (28) and (29) on the 
relative time s° = x°—y° is assumed to be determined by these quasiparticles. 
According to (39), given later, the Hartree-Fock energy is given by u p = 
£ + 2g[n c (x) + n(x)] (see also 0). 
Inserting the above approximations in (25) gives 

id x o + - g[n c (x) + 2n(x)]j ${x) = 2ig®(x) dz° J d 3 z 

3 

n 



X 



' d 3 pi 



i=i 



J(mv c -p 1 -p2+P3)(z-x)-i(e c -Lxj 1 -u>2+uJ3)(z -x °) 



(2tt) 3 , 

x[/ 2 (l + /s)/i-(l + / 2 )/s(l + /0], (30) 
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where /j means f(pi,X). In writing (30) the expansion (following [2]) 
$(z) ~ <|)^ 2 .) e "TM; c (2;-a;)-je c (2 -x ) - g appjjg^^ w ith v c and e c being the conden- 
sates atoms local velocity and energy, respectively. 

In taking the limit t° — > — oo one may approximate (again following the 
time integral in (30) by 



dz°e^-^ + ^"-^f ( ) ~ f(x°)n5(e c -u Jl -u 2 + u J3 ). (31) 



The space dependence of the distribution functions may be simplified using 
a gradient expansion around x giving f(X) ~ f{x). The i?-integral in (30) 
can then easily be performed and one finally obtains 



id x o + g[n c (x) + 2n(x)\ ) <r>(.r) 



ig 2 §(x) 



J d 3 p 1 d 3 p 2 d 3 p 3 5(mv c - pi - p 2 + Pz)5{tc -Ui-u 2 + lo s ) 



(27T) 

x[/ a (l + /s)/i-(l + / a )/ 3 (l + /i)]. (32) 

This equation is in agreement with the result derived in [2]. The physical 
interpretation is obvious: the first term in the integrand of the RHS corre- 
sponds to the gain term, i.e. it describes the collisions of two non-condensate 
atoms scattering into a non-condensate and a condensate state. Similarly, 
the second term describes the inverse process leading to a loss of conden- 
sate states. The collision term is shown pictorially in figure 2. Dashed lines 
correspond to the condensate and full ones to the non-condensate field. 





gain loss 

Figure 2: pictorial representation of the collision integral 



3 Propagator equation 

In the previous section we derived a kinetic equation for the condensate field 
$ in terms of the non-condensate distribution function /. In order to obtain 
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a complete set of equations we thus have to find the Boltzmann equation 
fulfilled by /. Our starting point is the propagator equation of motion derived 
within the -^-expansion of the 2PI effective action to NLO (see (36) in [T5] ) 



□ + m 2 + — [$ c (x)$ c (x) + Gcc(x, x)] j Gui.U. ij) iS a i,S c (.r - //) 



+<S> a (x)K b (x,y) -i j dz° j d s zD(x, z)G ac (x, z)G db (z,y), (33) 
with (see (34) in |Ej) 

D(x, y) = i^5c(x ~y) + 3^-^0(2/' x )®a(x) 

-i-^ J dz° J d 3 zG ab (x,z)G ab (x,z)D(z,y). (34) 

We now sum equation (33) over suitable values of the indices (a,b), in other 
words we take (1,1) + (2,2) + z(l,2) — z(2, 1). Next we perform the non- 
relativistic limit and the field redefinition as before, leading to 

A % 
{-id x o + — - gn{x)) < T^*(x)^(y) >= -5 c (x - y) + $*(x)K(x, y) 

-~J dz° J d 3 zD(x,z)[< T^*(x)^*(z) ><T^(z)^{y) > 

+ < T^*(x)^(z) >< T^*(z)ty(y) >] , (35) 

with 

n(x) = n c (x) + n(x), (36) 

D(x, y) = 2ig6 c {x - y) + 4g®*(x)K(y, x) + 4g®(x)K*(y, x) 
+2ig J dz° J d 3 zU(x,z)D(z,y), (37) 

and K given in (21). Note that the correct non-relativistic replacement of 
the Klein-Gordon operator is here 

(□ + m 2 ) — ► -2m(-id x0 + —), (38) 

2m 

as it acts on iff*. 
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3.1 Boltzmann equation 

One may now derive a Boltzmann equation for the non-condensate distri- 
bution function /. Starting from equation (35) we apply again the Popov 
approximation and we neglect contributions in higher order than g 2 , giving 

(~id x0 + - 2gn(x)) < T**(x)*(y) >= 

l -5 c {x-y)-2ig 2 dz° d 3 z (39) 



c 



2 

x [$*(a;)$(jz) < T^*(x)^(z) >< T^(x)^*(z) >< T^*(z)^(y) > 
+<F(z)$(x) < TV(z)V*(x) >< TV*(x)V(z) >< T^*(z)^(y) > 
+&*(x)$(z) < TV*(z)V(x) >< TV*(x)V(z) >< T^*(z)^(y) > 
+ < Tty*(x)^(z) >< T^(x)y*(z) >< T^*(x)ty(z) >< T^*(z)^(y) >]. 

Following, e.g., Blaizot and Iancu (see section 2.3 of [22]) and using the tech- 
niques described in the previous section we introduce the relative and center 
of mass coordinates s and X like before and perform a Fourier transform 
according to (28) and (29). 

As a result the Boltzmann equation finally becomes (with T = X°) 

[dr + Ijf • % - 2g(d s n(X)) ■ dp] f(p, X, T) = C 12 + C 22 . (40) 

The r.h.s. of the above equation consists of two collision integrals describ- 
ing two body collisions between non-condensate atoms (C 22 ) and collisions 
involving one condensate atom {C\ 2 ). They explicitly read 

2g 2 f 

C22 = 777-77 / d 3 p 1 d 3 p 2 d 3 p 3 5(p 1 + p 2 -p 3 - p)5{uji + lo 2 -u 3 -u) 
(2tt) 5 J 

x [hh{l + / 3 )(1 + /) - // 8 (1 + / 2 )(1 + fx)] , (41) 
4g 2 n c (X) f 

Ci 2 = 7—^ J d 3 p 1 d 3 p 2 5(mv c + pi - p 2 - p)S(e c + uj x - lo 2 - u) 
x[/ 1 (l + / a )(l + /)-// a (H-/ 1 )] 
J d 3 p 1 d 3 p 2 5(mv c - p[ - p 2 + p)5(e c - uj x - uj 2 + u) 
x[/i/ 2 (l + /)-/(l + / 2 )(l + /i)]. (42) 



2g 2 n c ( X) 

(27T) 



This result is in agreement with the one derived by Griffin, Nikuni and 
Zaremba |2J and by Stoof [Hj. 
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The approach followed in this paper provides a systematic approximation 
scheme for the non-equilibrium evolution of a Bose-Einstein condensate. 
Starting from the nonperturbative equations (25,35) for scalar fields including 
scattering and memory effects it gives as approximation a practicable system 
of kinetic equations (32,40) at order g 2 describing the dynamical properties 
of condensate formation. These equations have already been tested numeri- 
cally using N-body simulations [22] . 

A full solution of the coupled 2PI equations (25,35), however, requires exten- 
sive numerical work which is beyond the scope of this paper. 
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